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TECHNICAL NOTE 2282 

AN IMPROVED APPROXIMATE METHOD FOR CALCULATING 
LIFT DISTRIBUTIONS DUE TO TWIST 
By James C. Sivells 

SUMMARY 


A new method is presented for calculating the lift distribution due 
to twist which gives a much closer approximation than the empirical 
method of Schrenk (NACA TM 9^8) and requires about the same amount of 
computing. The new method, based on lifting-line theory, makes use of 
the lift distribution due to angle of attack and takes into account the 
aspect ratio of the wing. The twist may be that of the washout incor- 
porated in the wing design, of aeroelastic deformations, of deflected 
flaps or ailerons, or of downwash induced by another lifting surface or 
by the jet boundary of a wind tunnel. Examples are presented for the 
four main types of twist: symmetrical, anti symmetrical, continuous, and 

discontinuous. The applicability of the method to swept wih'g^ i's also 
presented. 


INTRODUCTION 


For many purposes it is desirable to calculate quickly an approxi- 
mate spanwise lift distribution on a wing due to its twist. The twist 
may be that of the washout incorporated in the wing design, of aero- 
elastic deformations, of deflected flaps or ailerons, or of downwash 
induced by another lifting surface or by the jet boundary of a wind 
tunnel. One approximation, which has been used extensively in the past, 
is that of Schrenk (reference l) modified by fairing out discontinuities 
as suggested in references 2 and 3. Although Schrenk’ s approximation may 
be adequate for some purposes, it is not as accurate as often desired 
since it does not take into account the aspect ratio of the wing and 
weights the twist angle arbitrarily according to the spanwise chord dis- 
tribution. The ability to obtain a closer approximation with about the 
same amount of computing is therefore desirable. The method presented 
herein has been developed to provide such an approximation. 

Gdaliahu, using lifting-line theory, has rigorously proved in refer- 
ence A that the average angle of attack of a twisted wing can be obtained 
by a spanwise integration of the local angle of attack weighted according 
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to the additional spanwise load distribution. The average angle multi- 
plied by the three-dimensional lift-curve slope gives the lift coeffi- 
cient. These results are proved herein in a somewhat different manner. 
It is also shown herein that the moment of the lift distribution of a 
twisted wing can be closely approximated by a simple expression. It is 
then reasoned that, if both the lift and moment given by the approxi- 
mation agree with the lift and moment obtained from more exact calcula- 
tions, then the approximate spanwise load distribution must agree fairly 
closely with the theoretical spanwise load distribution. A few examples 
given herein show this reasoning to be correct. 


SYMBOLS 


A 

A n 

Bn 


^n 

E 


E' 


a 


aspect ratio 

coefficients of trigonometric 

coefficients of trigonometric 
per unit angle of attack 

increment in lift coefficient 

coefficients of trigonometric 
due to twist 


series for lift distribution 
series for lift distribution 

due to twist 

series for lift distribution 


edge-velocity fact or for symmet rical lift 
distributions (\|l + il_co|£A_ > \ 


A2 / 


edge-velocity fact or for anti sym metrical lift 

l6 cos 2 A \ 


distributions 


(\F 


A c 


three-dimensional lift -curve slope per degree 


®o 

b 


two-dimensional lift-curve slope per degree 


section lift-curve slope per degree 



span of wing 


c local chord of wing 

c mean geometric chord (b./A) 
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c l 



c *b 

c le 

pb/2V 

y 

y* 

a 

°a 

e 

7 

6 * 

€ C 


e 

A 


section lift coefficient 

section lift coefficient for additional lift distribution 

section lift coefficient for additional lift distribution at 
a wing lift coefficient of unity- 

section lift coefficient for basic lift distribution 

section lift coefficient for lift distribution due to twist 

wing-tip helix angle, radians 

spanwise coordinate 

spanwise coordinate at position of discontinuity in twist 

angle of attack, degrees 

induced angle of attack, degrees 

angle of twist, positive for washin, degrees 

average angle of twist, degrees 

magnitude of discontinuity in twist, degrees 

correction factor used in fairing discontinuity in twist, 
degrees 

spanwise coordinate ^cos" 1 

sweep angle of wing quarter-chord line, degrees 
DEVELOPMENT OF METHOD 


The approximate method developed herein for calculating lift dis- 
tributions due to twist is based upon the determination of exact or 
approximate expressions for the area and moment of the lift distributions 
due to the twist in terms of the twist itself, the average twist, the 
additional lift distribution due to angle of attack, and the wing lift- 
curve slope. The lift distribution due to twist is designated as c^cyb 

which becomes equal to the basic lift distribution c Zfcc/ b when the wing 
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lift coefficient is equal to zero. In general. 


C Z € C 


Ci c c, c 


=A C L ^ + ^ 


where 


n 


ac l =. 


A 


L ¥ if) 


Reference 4 and appendix A prove that 



so that 



and the average twist 




( 1 ) 


(2) 


(3) 




(5) 


SYMMETRICAL TWIST DISTRIBUTIONS 

For symmetrical twist . distributions, the integral of equation (l) 
can he written as 
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by using the preceding relationships. Since the last integral of equa- 
tion (6) is the integral of the basic lift distribution, it is equal to 
zero and may be divided by any arbitrary constant without changing its 
value . Therefore , 



From a consideration of the moment of the lift distribution due to twist, 
appendix A shows that 



where 


F = 


1 + 


3&> a 
*2 A 


Since the areas are equal (equation (7)) and the moments of the same 
areas are approximately equal (equation (8)), it is reasoned that the 
distributions are approximately equal. Therefore, 


c l e c 
b 



€ - ? \ c *al c 

Fib 


a 


(9) 


A comparison of equations (l) acid (9) shows that the additional part of 
the lift distribution due to twist is given by the exact expression 



£ 


c *al c 

b 


a 


(id) 


and the basic part of the lift distribution due to twist is given by the 
approximate expression 




b 


'\> 


a 
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Since 


F = 1 + 


360 


a 


,2 A 


equation (ll) may be written as 


c 2-k c 


C7 c 
e - e *al 


^ _ 

b ^ + 360 a b 

+ rc2 A 

By use of the further approximation (see reference 5) that 


( 12 ) 


~ ~ A n ~ A jf2 

d -"v dm ~ —————— 

AE+2 AE+290 

equation (12) may be written as 

c *b c A(6 - e) c *al c 
~ ~ AE + 6 “V - a “ 


(13) 


Anti symmetrical Twist Distributions 

For anti symmetrical twists, AC^ and e are both equal to zero, so 
c 7 c c 7 c 

that = . From a consideration of the lift for one semi span of 

b b 

the wing, appendix A shows that 



where 


F' 


1 + 


180 a 
*2 A 
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From a consideration of the moment of the lift distribution, it is also 
shown in appendix A that 



Again, since the areas are approximately equal (equation (l4)) and the 
moments of the same areas are also approximately equal (equation (15)), 
it is reasoned that the distributions are approximately equal. Therefore, 



(16) 


or 



€ 



1§0 

2 


b 


a 


(17) 


which may be further approximated as 


c-, c c 7 c 

L e ^ Ac ‘■al 

b ~ AE + if- b 


( 18 ) 


The use of the edge-velocity factor E for symmetrical lift distributions 
is inherent in the method developed herein. For anti symmetrical twists, 
the use of the edge-velocity factor E* has been shown to give results 
which are in closer agreement with those obtained by lifting-surface 
theory. (See reference 5.) In order to take this effect into account, 
equations ( 17 ) and (l8) may be modified by multiplying the right-hand sides 

AE + if. AE + 2 

by the factor , or a fictitious lift-curve slope a’ = a — 

AE 1 + if- . AE' + 2 

may be used in equation ( 17 ) instead of the lift-curve slope a. 
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ILLUSTRATIVE EXAMPLES 


In order to illustrate the degree with which the lift distribution 
due to twist may be approximated by the method described, results of 
calculations are presented in figures 1 to 7 for the four main types of 
twist: symmetrical, anti symmetrical, continuous, and discontinuous. 

The wing chosen for four of the examples had a taper ratio of 0.5, an 
aspect ratio of 6.7k, and a corresponding edge-velocity factor of 1.043. 

To show the effect of aspect ratio, an additional example is presented 
(fig. 3) for a wing with a taper ratio of 0.5, an aspect ratio of 13.92, 
and a corresponding edge-velocity factor of 1.010. A two-dimensional 
lift-curve slope of 0.1097 was assumed. The additional .lift distribution 
(due to angle of attack) was obtained by the method of reference 6. Equa- 
tions’ (12) and (17) were used inasmuch as they yielded slightly closer 
approximations than equations (13) and (18) for the comparison with result 
obtained by the more nearly exact methods of references 6 and 7 utilizing 
linear section lift curves. Also presented are results obtained by 
Schrenk's approximation as described in references 1 to 3 hut modified 

by using the section lift-curve slope a 0 = -=r instead of the two- 

dimensional slope aoo. Therefore, 



(19) 


Symmetrical Twist Distributions 

Continuous twist .- The lift distributions shown in figure 1 are for 
the induced twist due to the jet boundary of a wind tunnel. The resulting 
distributions contain both the additional and basic parts. The main dis- 
crepancy between the two curves occurs at the root where the approximate 
method is unable to show the zero slope required by potential theory. 

The basic parts of these distributions are shown in figure 2 along with 
that calculated by Schrenk's method. Both approximate methods show the 
same discrepancy at the root but Schrenk's method shows a much larger 
discrepancy near the tip. The differences between Schrenk*s method and 
the "method of this paper are more clearly shown in- figure 3 for the wing 
of aspect ratio 13*92 and the same twist as for figure 2. 

Discontinuous twj st . - The lift distributions due to the deflection 
of 50 -per cent -span flaps are shown in figure 4. The direct application 
of the method will give a discontinuity in the lift distribution which 
is impossible according to potential theory. In this figure, the dis- 
continuity in twist distribution was faired with ellipses as indicated 
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or in terras of a correction factor to be added to the discontinuous 
twist distributions: 


For 


2 y* 


°< 17 < „ 

b b 



( 20 a) 


( 20 b) 


where e* is the magnitude of the discontinuity in twist at the span- 
wise position '2y*/b of the discontinuity. This type of fairing is 
arbitrary and some other type might be used to give better results. The 
discrepancies shown in this figure are wholly due to discrepancies in 
basic lift distribution as shown in figure 5* The unfaired distributions 
are shown for comparison. The root bending moment of the unfaired curve 
obtained by the method described herein is less than 1 percent different 
from that of the more nearly exact curve. Furthermore, if the ordinates 
of the faired curve are multiplied by a constant to give the same root 
bending moment as that of the unfaired curve, a much closer approximation 
will be obtained. The closeness of the approximation obtained by Schrenk’s 
method obviously depends upon proper fairing not only at the discontinuity 
but also near the tip. 


Anti symmetrical Twist Distributions 

Continuous twist . - The lift distributions of a wing rolling with 
a tip helix angle of 0.01 radian or O. 573 0 are shown in figure 6 . In 
this case the agreement between the * approximate' method and the more 
nearly exact method is extremely good, whereas Schrenk's method is not 
nearly so good. 

Discontinuous twist .- The lift distribution due to the deflection 

of 50-percent- span ailerons is shown in figure 7. As in the case of 
the partial-span flaps, both the unfaired and faired distributions are 
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shown) the root tending moment of the unfaired curve is practically 
identical with that of the more nearly exact curve. The discontinuity 
in twist was faired in the same manner as for the case of the 50 -percent - 
span flaps, the difference being that the signs of the correction factors 
were reversed. Again, if the ordinates of the faired curve are multi- 
plied by a constant to give the same root bending moment as that of the 
unfaired curve, a much closer approximation will be obtained. The inade- 
quacy of Schrehk's method is clearly shown. 


DISCUSSION • 
Comparison of Methods 


The examples shown were chosen to indicate some of the similarities 
as well as the differences between Schrehk's method and the method devel- 


oped herein. In the case of an elliptical wing for which 


C ^aI C _ c 




is readily evident that Schrehk's original approximation, which uses the 
two-dimensional lift -curve slope 



(V 


e-r-g c 

2 b 


a co 


will give the same result as the improved approximation for an aspect 

ratio of approximately 6 for symmetrical twists and for an aspect ratio 

of approximately b for anti symmetrical twists. An apparent improvement 

in Schrenk's approximation would be the use of the appropriate aero- 

A A 

dynamic induction factor, for symmetrical twists and 

AE + 6 AE* + 4 

for anti symmetrical twists, rather than the factor of Such a method 

would then become that of the so-called strip theory modified by the 
aerodynamic induction factor. The factor — - — however, should be 

used to obtain the basic lift distribution due to symmetrical twists 

rather than the factor — — commonly used. 

AE + 2 


In addition to taking into account the aspect ratio of the wing, the 
improved method makes use of the additional lift distribution instead of 
the chord distribution used by Schrerik. As the departure from elliptical 
lift and chord distributions becomes greater, the method developed herein 
becomes better as compared with Schrenk’s method. Although the additional 
lift distributions are not quite as readily available as chord distributions, 
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those of a large variety of plan forms can be found in reference 8. 
Even using the approximations of reference 1 or 9 for the additional 
lift distributions is better than using the chord distribution. Once 
the additional lift distribution has been obtained, the time required 
for computing the basic lift distribution is essentially the same for 
the improved method as for Schrenk's method.. 


Application to Swept Wings 

Inasmuch as the approximate method presented herein makes use of 
the additional lift distribution and three-dimensional lift-curve slope, 
it might be expected to yield satisfactory results for swept wings even 
though the method is based upon lifting-line theory. The adaptation of 
the formulas for swept wings is given in appendix B. Comparisons for 
two sweptback wings are presented in figures 8 to 10. For figure 8, the 
sweepback was , the aspect ratio was 1.5, and the taper ratio was 0.5* 
For figures 9 and 10, the sweepback was 60°, the aspect ratio was 3*5, 
and the taper ratio was 0.5* The additional lift distributions and 
three-dimensional lift-curve slopes were obtained from reference 10 for 
use in equations (12) and (17). Curves from references 11 and 12 are 
used for comparison. In the application of Schrenk’s approximation to 
the case of the sweptback wings, it was evident that using the two- 
dimensional II ft -curve slope would give exaggerated values for the lift 

/ \ a co cos A 

distribution. Equation (19) was therefore used, for which a q - » 

E 

in order to provide a more nearly equal basis for comparison. 

Figures 8 and 9 for linear symmetrical twist distributions indicate 
that the improved approximation does yield reasonably satisfactory results 
even for the extreme configurations selected. The fairly good agreement 
obtained in figure 9 with the modified Schrenk's approximation is for- 
tuitous in the same manner as in figure 2 since — was of the order 

cos A 

of 6 for this configuration. For the configuration of figure 8, the value 

of - - was much less, and Schrenk's approximation is in error in the 
cos A 

opposite direction from that shown in figure 3 for a high aspect ratio. 

Figure 10 for the anti symmetrical twist distribution of a rolling 
wing indicates less satisfactory results. The general shape and lateral 
center of pressure are in agreement with reference 12, but the magnitude 
is about 18 percent too great. Equation (17) was used directly for this 
calculation. The suggested modification of using the fictitious lift- 
curve slope a' would decrease the error to about 10 percent. This error 
may be too large to be tolerated for the computation of the damping-in- 
roll derivative but may be sufficiently close for other purposes in which 
anti symmetrical twists are involved. 
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CONCLUDING REMARKS 


An improved approximate method has been developed for calculating 
the basic lift distribution due to any of the four main types of twist: 
symmetrical, antisymmetric al, continuous, or discontinuous. The method 
makes use of the additional lift distribution due to angle of attack 
and takes into account the aspect ratio of the wing. Although the 
method is based upon lifting-line theory, it gives reasonably good 
results even for swept wings. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., May 2k y 1950 
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APPENDIX A 

DERIVATION OF APPROXIMATIONS FOR AREA AND MOMENT OF 
LIFT DISTRIBUTIONS DUE TO TWIST 
General Equations from Lifting-Line Theory 


In the following derivation, lifting-line theory is used and the 
section lift curves are assumed to be linear. The spanwise lift dis- 
tribution is expressed as the trigonometric series (reference 6 ). 


c l c 

b 


" 2Z 


sin n 0 


(Al) 


where 6 is defined by the . relation cos 6 = The lift distribution 

is also related to the section lift-curve slope by the expression 


c z c 

b 


a o c / >. 

— (fL - ai) 


(A2) 


where the induced angle of attack, in degrees, is expressed by 


aq = — V"" nA n sin nQ (A 3 ) 

kit sin Q c — 

The lift distribution can be divided into two distributions by the 
substitution of A n = aBn + C n in equation (Al) . This substitution 
gives 

= a ®n sin n0 + C n sin n0 (hk) 


where ^ sin n 0 is the lift distribution per unit angle of attack 
and ^n sin rid is the lift distribution due to twist. The additional 
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lift distribution (that due to angle of attack for a wing lift coefficient 
of unity) is 



— V Bjj sin n 6 
a — 


The lift distribution due to twist is designated 


(A5) 


c l c 

— 1_ = c n sin n0 

"K t 


An integration of equation (A5) along the span yields 


/-I 


c 7 c 
z al 


I n 

y~ B n sin n0 sin Q d0 


or 


p 7tB- 

A “ 2a" 


■R 

Bi = — 


The increment in lift coefficient due to twist is 


ac l = | 


c 7 C 


-1 


(?) 


r 


C n sin n0 sin 6 d 6 


JO 


(A 6) 


(A7) 


(A8) 


ac l 


AjtC]_ 


k 


(A9) 
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Substituting equation (A8) into equation (A9) yields 


ACl = a 


Bl 


The average twist is 


f _ ac l _ Cl 

■ ~ - 51 


For a wing with twist, equation (A2) may be written 


c z c _ a 0 c 


(a + e - aj 


b b 

which can be combined with equations (Al) and (A3) to give 


XAn 


sin n 0 = ^2. ( a + € - T — 

b V 4 tc sin 0 


y~ nA n sin n©j 


Equation (A13) can be divided into the two equations: 


Y~ B n sin n0 = -2— / 1 


r n Bn sin n0) 

sin G J 


and 


F C m sin m0 = —2— ( e - — — mC^ sin m0 ) 

b ^ n: sin 0 y 


where n has been replaced by m for subsequent use. 
can be rearranged in the following manner: 


Equation 


e sin 0 - — Y" mC™ sin m0 
jt • 


(A10) 

(All) 

(A12) 

(A13) 

(A1»0 

(A15) 

(A15) 


~r Cm sin m0 sin 0 


(A16) 
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Multiplying equation (Al4) by equation (Al6) yields after rearranging: 


€ ^ Bn sin n0 sin 6 ~ y C m sin m0 sin 6 + 


(m - n-)B n C m sin nQ sin m0 (AIT) 

71 n m ■ 

which is the basic equation used hereinafter. Only odd values' of the 
index n exist because the lift distribution given by equation (A5) is 
symmetrical. Values of the index m are odd for lift .distributions due 
to symmetrical twist distributions and even for lift distributions due 
to antisymmetrical twist distributions. 


Symmetrical Twist Distributions 

For symmetrical twist distributions, the same result will be obtained 
whether equation (A17) is integrated along the entire span or along the 
semispan. Integrating the double-summation term of equation (A17) gives 



^5 y 
* n 


E (m - n)B n C m sin n0 sin m6 d0 = 0 
m 


so that 

JO 


t/2 

6 E Bn sin nS sin 8 dS = I E C m sin m0' sin 8 dS (Al8) 



(A19) 


which was proved in a different manner in reference 4. Equation (A19) can 
be readily shown to be valid for asymmetrical twist distributions if the 
integrations are performed along the entire span. If € is replaced 
by e in equation (Al8), the integration yields 


TT El = £ Cl. 


_ Jt 


which is consistent with equation (All). 



Subtracting £ sin sin ® from both sides of equation (A17) gives 
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ci 

o 0) 

•p P rd 
d p 
0 ) 0 ) 
d S 


c! 

O <D 
•H ^ O) 

p p r d 

p bOQ> 
•H d 
U *H W 
P > ; O 


•H Pi 'd 
H -H <U 

HI M H 
43 *3 Pf 
P O -H 
•H P 
Pi P H 
odd 

1 — I H 
P 0) 

p c< w 

CO -H 
ft aJ 


lT\| 1— t 
pq I CVJ 


col _ 
pq I ir\ 


1 — I rH 

o pq 


a <y 
o > 

•M *H 
P W> CM 
«J 

d w k 
o« P C_ 

QJ O 


■“ Is .CVI 


H H 
O W 


n)BnC m sin n 0 sin s 6 co« 0 d0 = ^ £^03 - B 3 c l) 2r( B l c 5 “ B 5 c l 
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where 


F « 


1 + 


90 


Bl 


= 1 


360 a 
2 A 

IX 


Anti symmetrical Twist Distributions 

For anti symmetrical twist distributions, the average twist is 

ci c c^, c 

equal to zero, so that g = For this case, the integration of 

equation (A17) along the entire span is zero. In order to obtain a 
relationship involving the lift of one semispan, equation (A17.) is 
integrated along the semispan. Thus, 


Prt/2 


e ZZ ^ sin n 0 sin 0 d0 = 


10 


'it/ 2 


C m sin m0 sin 0 d0 + 



it 


5T H (m - n)BnC m sin n0 sin m0 d0 
n m 


(A24) 


The integrations of the right-hand terms of equation (A24) yield 


Vr/2 


YZ sin sin Q = 


2C 2 


4C4 + 6C5 + 

15 
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and 


Vc/2 


— XZ XI ( m - n )®n < -'in s ^- n s i n ^ = 

* n m 




2B 3 2B 5 

+ — + — + ... 


it \ 5’ 7 9 


*>C 6 t± + 6 -h + 6 »l + ...) + ... 

« b \7 9 11 / 


Therefore, 

rv/2 


/ 2Co ^k 6c 6 W kR 

e n B n sin n0 sin 0 d0 = [ — - + — + . .. f(l + — B l) + 


Jo 


45 

jr 


C 2 


/2B 3 2B. \ /8B-L 

W’t 5 * -J-Hr 


3 15 35 


4Bo ^-Btr 

+ + + 

7 9 


••■) 


06 


(2b^ + 6B 1+ 6^ 


35 9 . 11 


■) 


+ . . . + 


(A25) 


If the last term (in brackets) of equation (A25) is neglected, equa- 
tion (A24) becomes 

'jt/2 P*/2 

B n ' sin n0 sin 0 d0 % j XI 0 m sin sin ® 


JO 


or 


'1 f’l 

CtC I Ct C 

e l al 


- l * d /2y 


Jo 


b 


(?) 


F* b 


*•(?) 


(A26) 
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where 


F‘ = 1 + 22 Bi 

it 


= 1 + 


is 

n 

l80 a 
A 


In order to obtain a relationship involving the moment of the basic 
lift distribution, each term of equation (A17) is multiplied by cos Q d0 
and then integrated along the semispan. These operations give 


n */ 2 Pit/ 2 

e Y1 B n sin sin 0 cos 0 d0 = | y~ Cm sin m0 sin 6 cos 9 d0 + 


00 


V/2 


(m - n)B n C m sin n0 sin m0 cos 9 d0 


do 


n m 


(A27) 


= 8 ° 2 + t[( Bi ' B 3) C 2 + ( b 3 ' B 5) c b + (b 5 - b 7 )c 6 + ...] 


Therefore, 


V/2 


e E ^ sin n0 sin 0 cos 0 d0 = 


00 


| c 2 (l + ^ Bj) - |[b 3 (c 2 - c u ) + B 5 ( Cl , - Cg) + ...J (A28) 

If the last term (in brackets) of equation (A28) is neglected, equa- 
tion (A27) becomes 


fV/2 P«/2 

— YZ ^n sin sin ® cos &d9z C m sin m0 sin 0 cos 0 d0 

00 F ‘ Jo 
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APPENDIX B 

ADAPTATION OF FORMULAS FOR SWEPT WINGS 


For swept wings, both the edge -velocity factors and the three- 
dimensional lift -curve slope must be modified to include the effects of 
sweep. The edge -velocity factors are listed in the section entitled 
"Symbols." The formulas as listed are different from those given in 
reference 5 in that the cos 2 A term used herein is given as a cos^A 
term in reference 5* Although there Is no rigorous proof that either 
factor is correct, the use of the cos^A term has been found to yield 
answers in better agreement with those given in references 10 and 12. 
For example, the 10-percent error mentioned previously in connection 
with figure 10 would be increased to 1 6 percent by using the cos^A 
term. The three-dimensional lift-curve slope can be written 


tions can be written 


“ + 2 
cos A 

*00 

A 

*2 

AE +2 
cos A 

90 

slope used 

for a 

AE 

+ 2 

cos A 

a AE’ 

+ p 

cos A 


.... A 


AE ' , 

cos A 

d 00 

2 

A 



AE 


+ 2 


90 


cos A 
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As a matter of interest, using the edge-velocity factors as defined 

herein ‘yields the result that the term — approaches the value of 2 

cos A 


AE 1 

and the term approaches the value of 4 as the aspect ratio 

cos A 

approaches 0. This result is in accordance with lov-aspect-ratio theory 
(reference 13 ) • ■ 

Using the relationships for a and "a' in equations (12) and (17) 
yields the following equations: 

For the symmetrical twists, 


or 


C H) C ^ g - € C *al C 
b * 36 O a b 

tc 2 " A 


(Bl) 


C V 

A( € € ) c 2al c 

b * 

AE , ^ b 


cos A 

where 


r = a 

I • <*) 

For anti symmetrical twists, 


c 7 c 

.. Z b _ 

_ c 7 c 

1 € l al 

b ~ ■ 

l 80 a 1 b 
«2 A 

or 


c l c 

As *al 

-rt-vr; a. 


cos A 


(B2) 


(B3) 


(B4) 


(B5) 


For discontinuous twists, the discontinuity must be faired out. For 
convenience, this fairing may be accomplished through the use of ellipses 
as indicated in equations (20). 
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Figure 6.- Basic lift distribution for anti symmetrical continuous twist 
due to rolling. — = o.Ol radian; taper ratio, 0.5; aspect ratio, 6.7^; 
edge -velocity factor, 1.0^3. 















Figure 8.-. Basic lift distribution per degree of linear twist. Sweep, k-5' 
taper ratio, 0.5; aspect ratio, 1.5 • 













Schrenk's method (modified) 



Figure 9*- Basic lift distribution per degree of linear twist. Sweep, 60° 
taper ratio, 0.5; aspect ratio, 3 . 5 . 
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